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Recent experimental and theoretical investigations on asymmetric field emission induced by weak 
femtosecond laser pulses and also its emission mechanisms are briefly reviewed. The emission mech- 
anisms are discussed further for a wider range of DC fields and laser power. It appears that firstly 
photo-assisted field emission from lower-excitation order grows in the higher DC fields and secondly 
our simulations can be applied only for lower laser power. 



I. INTRODUCTION 

Field emission is the electron emission from a solid to 
the vacuum via electron tunneling through the poten- 
tial barrier at the surface. It is driven by applying a 
strong electric field (2-5 V/nm) to the solid [IQ. Field 
emission from metallic tips with nanometer sharpness has 
been particularly interesting because its electron beam 
is highly coherent and bright due to a nanometer- or 
atomically- sized emission area, and it has been widely 
used in both basic research and applications 

Field emission from such a tip can be induced with 
a moderate DC field and laser irradiation onto the tip 
apex. So called laser-induced field emission was intro- 
duced some time ago, first by using a continuous laser 
beam Q and later by pulsed lasers [ltj E2 ■ More re- 
cently, ultrashort pulsed field emission by using femtosec- 
ond laser pulses was realized 

mm. The p uised fieid 

emission can be emitted from almost the same emission 
area as that of field emission [lj|- Hence potentially a 
spatio-temporal resolution down to the single atom and 
attosecond range appears to be possible @, H, EH Elf. 
Moreover, due to plasmonic effects, emission patterns 
were found to be asymmetric and can be controlled by 
laser polarization. In effect, an ultrafast pulsed field- 
emission source with site selectivity on the scale of a 
few tens of nanometers was realized [13, EH- So far, 
the emission mechanisms have been intensively investi- 
gated for weak and strong field regimes [12h27| . Such 
electron sources would be very attractive for both basic 
research and new applications like time-resolved electron 
microscopy, spectroscopy, holography, and also free elec- 
tron lasers pj, EI EI M, ^2MM ■ 

The laser-induced field emission is considered to be 
driven through the following two steps. 1. Illuminating 
the metallic tip with the laser pulses creates enhanced 
optical electric fields at the tip apex due to plasmonic 
effects (local field enhancement) [3lU33l |. 2. The en- 
hanced fields induce pulsed field emission in combination 
with a moderate DC voltage applied to the tip. Depend- 
ing on the strength of the enhanced fields, different field 
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emission mechanisms become dominant jl2H27| . For rel- 
atively weak fields, single-electron excitations by single- 
and multi-photon absorption are prevalent, and photo- 
excited electrons are tunneling through the surface po- 
tential barrier (photo-assisted field emission) or emitted 
over the barrier (photoemission) . On the other hand, 
very strong fields largely modify the tunneling barrier 
and prompt field emission from the Fermi level (optical 
field emission) . 

In this manuscript, our recent investigation on the op- 
tical control of field emission sites and its emission mech- 
anisms are briefly summarized in sections 2 and 3, respec- 
tively. Optical control of field emission sites was realized 
because the above mentioned local field distribution on 
the tip apex becomes asymmetric, thus driving the asym- 
metric field emission [bil . EH . Energy distribution curves 
(EDCs) of the emitted electrons revealed that their emis- 
sion mechanism is mainly dominated by photo-assisted 
field emission poj |. In section 3, we extend our discussion 
on EDCs for a wider range of DC fields and laser power. 
The results show a change of the emission processes and 
the spectral shapes with varying both parameters, and 
the validity of our simulations for EDCs appears to be 
limited to lower laser power. 



II. OPTICAL CONTROL OF FIELD EMISSION 
SITES 

In our previous work, we have investigated electron 
emission patterns induced by femtosecond laser pulses 
from a clean tungsten tip apex, and compared them with 
those of field-emitted electrons without laser excitation 
as schematically depicted in Fig. 1. A tungsten tip with 
its axis along the (Oil) crystal direction is mounted in- 
side a vacuum chamber. Laser pulses are generated in a 
Thsapphirc oscillator (wave length: 800 nm; repetition 
rate: 76 MHz; pulse width: 55 fs; average laser power 
Pl- 20 mW) and introduced into the vacuum chamber. 
The laser light was focused to 4 [im in diameter onto 
the tip apex, and emitted electrons were detected by a 
two-dimensional detector. The polarization vector can 
be changed within the transversal (x, z) plane by using 
a A/2 plate. The polarization angle 8p is defined by the 
angle between the tip axis and the polarization vector; 
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details are depicted elsewhere [13, [l5| . 

We observed a striking difference in symmetry of the 
two patterns. Without laser, a typical field emission 
pattern of a clean W tip was observed. The most in- 
tense electron emission was observed around the (310)- 
type facets and relatively week emission from (lll)-type 
facets as indicated in the picture. The intensity map 
roughly represents a work function map of the tip apex: 
the lower the work function is, the more electrons are 
emitted [HI, [34], HH . On the other hand, with laser ir- 
radiation, emission sites were the same as those without 
laser, but the emission pattern becomes strongly asym- 
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metric with respect to the shadow and exposed sides of 
the tip to the laser pulse. Furthermore, varying the laser 
polarization angle changes these distributions substan- 
tially as shown in Fig. 2(a). In effect, we have realized 
an ultrafast pulsed electron source with emission site se- 
lectivity of a few ten nanometers. Note that the emission 
site can be also selected by changing the azimuthal and 
the polar orientation of the tip apex relative to the laser 
incidence direction [3, [lj| . 

The physics behind driving asymmetric field emission 
is the asymmetric local field distribution on the tip apex. 
When a laser pulse illuminates the metallic tip, surface 
electro-magnetic (EM) waves like surface plasmon po- 
laritons are excited j36| , which propagate around the tip 
apex. Due to the resulting interference pattern, the elec- 
tric fields show an asymmetric distribution over the tip 
apex. To simulate the propagation of these surface EM 
waves and the resulting field distributions, we used the 
software package MaX-1 for solving Maxwell equations 
based on the Multiple Multipole Program [37H39| . A 
droplet-like shape was employed as a model tip as shown 
in Fig. 3(a), with a radius of curvature of the tip apex of 
100 nm, which is a typical value for a clean tungsten tip. 
The dielectric function of tungsten at 800 nm was used 

Fig. 3(a) shows the calculated time- averaged field dis- 
tribution over a cross section of the model tip, where the 
polarization vector has been chosen parallel to the tip 
axis (Op — 0°). The laser propagates from left to right. 
The field distribution is clearly asymmetric with respect 
to the tip axis, with a maximum on the shadow side of 
the tip. This is consistent with our observations in Fig. 1 
where the field emission is enhanced on the shadow side. 
Fig. 3(b) shows front views of time-averaged field distri- 



Fig 1: A schematic diagram of the laser-induced field emis- 
sion. The emission pattern is observed by a two-dimensional 
detector. The inset shows the field emission pattern without 
laser excitation. 
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Fig 2: (a) Observed and (b) simulated laser-induced field 
emission images for different light polarization angles Op. Op 
is defined by the angle between the tip axis and the polar- 
ization vector: Bp = is when the tip axis is parallel to the 
polarization vector. 



(b) 





8 P =0 



3 P =0 9 P =30 6 P =60 9 P = 90 6 P = 120 6 P =150 



Fig 3: The calculated time-averaged field distribution around 
the model tip is shown in a linear color scale in (a) for Op = 0° . 
Highest field values are represented in yellow. In (b) the time- 
averaged field distributions are given in a front view of the 
model tip for different polarization directions Op. The laser 
propagation direction is indicated by red arrows, (c) shows 
a schematic diagram of photo-assisted field emission from a 
nonequilibrium electron distribution.. 
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bution maps for different polarization angles. The field 
distribution changes strongly depending on the polariza- 
tion angle. 

From the calculated local fields, we further simulated 
the laser-induced field emission images by supposing the 
photo-assisted field emission mechanism as schematically 
depicted in Fig. 3(c). The current density of field emis- 
sion can be described in the Fowler-Nordheim theory 
based on the free-electron model 0, 0, SH S3- The 
simulated emission patterns in Fig. 2(b) are in excel- 
lent agreement with the observations shown in Fig. 2(a). 
This comparison clearly demonstrates that the observed 
strongly asymmetric features originate from the modula- 
tion of the local photo- fields at the tip apex [3 [ll| . 

III. MECHANISM OF LASER-INDUCED FIELD 
EMISSION 

A. EDCs of field emission and photo-assisted field 
emission 




Electron Energy (eV) Electron Energy (eVJ 



Fig 5: (a) EDC of field emission from the tungsten tip to- 
gether with a schematic diagram of field emission from a 
Fermi-Dirac distribution tunneling through a surface barrier. 
The transmission probability depends on the hatched area. 
The tip voltage Vu P was -2300 V, and the pinhole-plate volt- 
age Vpin was V. (b) shows experimental EDCs and schematic 
diagrams of photo-assisted field emission together with the 
magnified spectrum. The spectrum is taken at Vu p = -2300 
V, Vpin — V and the laser power Pl = 50 m W. 



In the previous section, the photo-assisted field emis- 
sion is shown to simulate well laser-induced field emission 
patterns. In order to confirm experimentally the emis- 
sion mechanism, we upgraded our previous field emission 
setup by implementing an electron energy analyzer (see 
Fig. 4) and measured EDCs of the emitted electrons. 
In order to distinguish electrons emitted from different 
facets of the tungsten tip, a plate with phosphor coating 
and pinhole in the center was installed as an aperture and 
counter electrode between the tip and the analyzer. The 
tip and the pinhole plate can be negatively and positively 
biased for field emission, respectively. The field emission 
pattern of the clean tungsten tip can be observed on the 
phosphor plate where the most intense electron emission 




Hemispherica 
Analyzer 



Fig 4: Schematic diagrams of the experimental setup. The 
emission pattern can be observed at the phosphor plate as in- 
dicated by green areas. An emission site is selected by the pin- 
hole (black dot) and EDCs are measured by a hemispherical 
analyzer. V P i„ and Vu p represents the pinhole-plate voltage 
and tip voltage, respectively. 



come from (310)-type facets. The pinhole is positioned 
at the edge of a (310) type facet as in Fig. 4. Note 
that the selected site is also the most intense emission 
site in the laser-induced field emission because linearly 
polarized laser light with the polarization vector parallel 
to the tip axis was used [f4|, [T^ . 

Fig. 5(a) shows an EDC of field emission at a tip volt- 
age Vup of -2300 V and a pinhole-plate voltage V P i n of 
V. The peak of the spectrum is defined as the Fermi 
energy Ep at eV [43|. The spectrum shows a typi- 



cal asymmetric peak, which can be understood by the 
diagram in Fig. 5(a). The field emission current is in- 
fluenced by two factors: 1) the electron occupation num- 
ber and 2) the transmission probability through the sur- 
face barrier [l], 0, SH S3- The occupation number is 
given by an electron distribution function f(E), which is 
the Fermi-Dirac distribution function in the case of field 
emission, and the transmission probability depends ex- 
ponentially on an area of the surface barrier above an 
energy of the emitted electron indicated by the hatched 
area. Therefore, the positive energy side of the spectrum 
falls off due to a rapid decrease of the occupation num- 
ber, while the negative energy side falls off because of the 
exponential decay of the transmission probability due to 
the increase of the surface barrier area. Thus a typical 
field emission spectrum shows such an asymmetric peak. 

Upon low-intensity laser irradiation, the electron dis- 
tribution is modified by single-electron excitations due to 
multi-photon absorption, resulting in a nonequilibrium 
distribution characterized by a step-like profile as illus- 
trated in the inset of Fig. 5(b) [H SI SI; the width of 
each step corresponds to the photon energy hu (= 1.55 
cV) . Here we identify the step edges of one-photon exci- 
tation (IPE) and 2PE as shown in the inset. In a real 
situation, the excited electrons relax mainly by electron- 
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electron (e-e) scattering on a time scale of a few femtosec- 
onds, and such an electron dynamics results in a smeared 
electron distribution as in the diagram. These features 
should be reflected in the EDCs. 

Fig. 5(b) shows an EDC of laser- induced field emis- 
sion at Vup = -2300 V and a laser power Pl of 50 mW. 
The spectrum shows the field emission peak undisturbed 
with identical shape and intensity as in Fig. 5(a), and 
additional peaks at the 1PE and 2PE edges are clearly 
observed. The latter show the same asymmetric shape 
as the field emission peak. Thus, photo-assisted field 
emission is confirmed experimentally. Regarding the rel- 
ative intensities, photo-assisted field emission from 2PE 
is much stronger than that from 1PE even though the 
occupation number at 1PE is higher; this is because the 
transmission probability at 2PE is quite high. Note that 
photo-assisted field emission from 2PE has not been ob- 
served for excitation with a continuous- wave laser Q: 
those electron distributions are supposed to be largely 
different from our case. 



with setting Vu p at -1500 V and varying V P i n from -1400 
V to 3000 V while Pl was fixed at 50 mW. These spectra 
show a smooth parametric transition between field emis- 
sion, photo-assisted field emission and photoemission. In 
the lower voltages (high negative), the tunneling proba- 
bility becomes suppressed and the emission spectrum is 
dominated by photoelectrons excited in a three-photon 
process as shown in the bottom spectrum of Fig. 6. The 
peak shows a spectral shape completely different from 
that of field emission. The peak maximum is located 
at approximately 0.65 eV below the 3PE edge. Since 
the transmission probability is unity throughout the pho- 
toemission regime, the peak shape reflects more closely 
the electron distribution function. Therefore the spectral 
shape indicates a strong modulation of the electron dis- 
tribution due to e-e scattering processes. For the higher 
voltages, the tunneling probability becomes larger and 
the photo-assisted field emission process becomes domi- 
nant. At around V p i n = 1000 V, photo-assisted field emis- 



sion from 2PE is dominant, and at higher V„in : 1 
assisted field emission from 1PE starts to grow [46 



ihoto- 



B. DC field and laser power dependences of EDCs 

The emission processes and the spectral shapes can be 
tuned by changing the DC fields. Fig. 6 shows various 
EDCs for different DC fields, where the EDCs are taken 
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Fig 6: Various EDCs for different V pin from -1400 V to 3000 
V. Vup is set at -1500 V. Magnification is also shown with 
each spectrum. 



The emission processes and the spectral shapes can 
also be tuned by changing the laser power. The laser 
power dependence of the EDCs was investigated by sys- 
tematically measuring EDCs for various tip voltages and 
laser powers as shown in Fig. 7. These spectra are nor- 
malized at their maximum intensities to focus on their 
spectral shapes. In this case, V p i n is grounded and only 
Vup is varied. As a rule of thumb, with increasing laser 
power, the electron emission from higher order excitation 
becomes dominant, and the spectra become broader and 
smeared out, which holds for each voltage setting. The 
reason why the spectra become smeared out with increas- 
ing laser power is that more excited electrons can scatter 
into lower energy states. 



The laser power dependence of intensities in the EDCs 
has also examined. Figs. 8(a)-8(d) show the raw data of 
Fig. 7 as solid circles in a semi-logarithmic plot for each 
voltage setting. The intensities of the photo-assisted field 
emission grow with increasing laser power for all voltage 
settings, while the field emission peak remains constant 
as shown in Fig. 8(a). From these raw spectra, inte- 
grated intensities for each excitation order are extracted 
as a function of laser power and shown as solid symbols 
in Figs. 8(e)-8(h). These intensities were obtained by 
integrating intensities over the respective energy regions: 
field emission (FE) is for -1.2 eV ~ 0.3 eV, photo-assisted 
field emission from 1PE for 0.3 eV - 1.85 eV and 2PE 
for 1.85 eV ~ 3.4 eV and photoemission from 3PE for 
3.4 eV ~ 4.95 eV. The solid lines are fitted curves with 
power functions. Exponents for each excitation order are 
also listed in the insets. These exponents are quite close 
to the correspondin g ex citation order, which is consistent 
with previous work [ill I20I? ]. 
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Fig 7: Measured and simulated EDCs for different tip voltages (DC fields) and laser powers (excitation constant C e xc', appears 
in eq. (f) of Ref. [2(|). The work function $ of 4.592 eV was used for the simulation. All the spectra are normalized at their 
maximum values. The measured EDCs are shown as green dots, and the simulated EDCs as red lines. 



C. Simulations of EDCs 

The experimental EDCs were reproduced and the in- 
fluence of electron dynamics on the EDCs was also ex- 
amined by simulating transient EDCs, taking into ac- 
count all relevant electron scattering processes, contain- 
ing both, contributions from new single-photon excita- 
tions from the current time step, and from the partly re- 
laxed electron distribution from all previous time steps. 
A laser pulse with a temporal width of 100 fs was moved 
across the emission site on the tip apex in 0.2 fs steps. 
The temporal evolution is covered from -200 fs to 1200 
fs for the laser power below 100 mW, while -400 fs to 
2400 fs for 200 mW and -400 fs to 3800 fs for 300 mW. 
The time zero is defined when the maximum of the laser 
pulse meets the emission site. At each time step, the 
transient electron distribution function and the result- 
ing EDC were calculated 0, H, SjSl. At the end 
of the each time window, the transient total yield be- 
comes approximately one percent of its maximum. All 
these transient EDCs were integrated over the complete 
time window. EDCs of field emission without the laser 
excitation was also calculated for the rest of one period 
of the laser pulse cycle (approximately 10 ns), and was 
added to the EDCs for the first laser-induced-emission 
period. The resulting time-integrated EDCs were nor- 
malized at the maximum intensity. Thus obtained sim- 
ulations were compared with the normalized experimen- 
tal EDCs in Fig. 7; comparison of the absolute intensi- 
ties were difficult because of the uncertain transmission 



of electron beams in the spectrometer. There are only 
three fitting parameters: 1) the work function $, 2) the 
DC field Fbci an d 3) the excitation constant C exc which 
is proportional to the laser power; the definition of C exc 
is given in the Ref. [2(| ■ More explanations on the sim- 
ulation is also provided in Ref. |2G |. 

The detailed fitting procedure is the following. Firstly, 
three fitting parameters are adjusted at one particular 
setting of the tip voltage and laser power (Vu p ~-500 V; 
Pl = 50 mW), and then time-integrated EDCs for the 
other settings were simulated by scaling up or down Fuc 
and C exc according to the corresponding tip voltage and 
laser power. This procedure was iterated until reason- 
able fitting was obtained for all settings. Thus all the 
experimental EDCs in Fig. 7 were reproduced by such 
simulations. Throughout all the various conditions, the 
simulations are in very good accordance with the spectral 
shapes of the experimental EDCs, as shown by solid lines 
in Fig. 7. Thus we conclude that the electron dynamics 
play a significant role in photo-assisted field emission. 

Also the intensities of the simulated EDCs were also 
compared with the experimental data, which arc shown 
as solid lines in Figs. 8(a)-(d). Because it is difficult 
to simulate absolute intensities, we renormalized all the 
simulated spectra in Fig. 7 at the maximum intensities 
of the raw experimental EDCs at 50 mW in each voltage 
setting. The integrated intensities for each excitation 
order are also extracted from those simulated spectra and 
shown as open symbols together with their fitting curves 
as dashed lines in Figs. 8(e)-(h). The simulations show 
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Fig 8: Semi logarithmic plots of raw experimental EDCs in Fig. 7: (a) V t4p =-2300 V, (b) V t i P =-1500 V, (c) V t i P =-1000 V, and 
(d) Vu p =-500 V. The solid lines indicate the simulations renormalized at maximum intensity in the EDCs at 50 mW of each 
voltage setting, (e) - (h) show the measured laser power dependence of intensities for each excitation order: (e) Vti p =-2300 
V, (f) Vti P =-1500 V, (g) Vii P =-1000 V, and (h) Vti P =-500 V. These intensities were obtained by integrating intensities over 
the respective energy regions: field emission (FE) is for -1.2 eV ~ 0.3 eV, photo-assisted field emission from 1PE for 0.3 eV 
~ 1.85 eV and 2PE for 1.85 eV ~ 3.4 eV and photoemission from 3PE for 3.4 eV ~ 4.95 eV. The solid symbols indicate the 
experimental data, while open symbols are from the simulations. The solid lines (experiment) and the dashed lines (simulations) 
are fitting curves with power functions; their exponents are also shown in the insets. 



quite good agreement with experimental data up to 100 
mW as can be especially seen in Fig. 8(a) and 8(e). 
The simulated exponents for each excitation order also 
show good accordance. Note that these exponents are not 
exactly 1, 2 and 3, but slightly higher. This is because the 
higher-order excited electrons decay due to the electron 
scatterings and some of the decayed electrons are emitted 
from the lower-order excitation region. 

Although simulations show good agreement up to 100 
mW, the simulations largely deviate from experimental 
data at 200 mW and 300 mW. These deviations become 
bigger with increasing laser power as can be especially 
seen in Figs. 8(d) and 8(h), and the simulated exponents 
become higher than the experimental ones due to the in- 
clusion of these data points. Because the experimentally 
obtained exponents are very close to the corresponding 
excitation order, we believe our simulation fails for the 
higher laser power regime. We speculate that the devi- 
ation is due to a lack of energy flow to the outside of 
the system (energy dissipation) in our simulation model. 
Without the energy dissipation, a decay speed of the ex- 
cited electrons would be slower than the actual decay 
[44]] . As a consequence, the electron emission lasts for 



a relatively long time, thus leading to an overcstimation 



T 




Time (fs) 

Fig 9: Temporal line profiles of the simulated electron pulses 
for Cexc = 1.7 x 10" 5 , 4.25 x 10" 5 , 8.5 x 10" 5 , 17 x 10~ 5 
and 25.5 x 10~ 5 . F DC is 0.502 V/nm for all the curves. The 
spectra are normalized at their maximum. 
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of the emission intensity. In our simulation, temporal 
line profiles of the simulated electron pulse as in Fig. 
9 show that the electron emission lasts much longer at 
Cexc = 17 xl0~ 5 , 25.5 x 10~ 5 than those at the lower C exc 
of approximately 200 fs. This indicates that the temper- 
ature of the phonon system rises significantly at higher 
laser power. According to the Ref. [HI, the inclusion 
of energy dissipation processes, especially ballistic trans- 
port, reduces the heating of the phonon system. Such an 
energy-dissipation term is necessary for a system where 
the finite optical penetration depth leads to a spatially 
inhomogeneous excitation 4J], which is also the case in 
the laser- illuminated tip apex (50j . Therefore, the inclu- 
sion of the energy dissipation in the simulation would 
solve this problem. Further investigation is required on 
this point. 



IV. SUMMARY 



with increasing the DC fields. For lower DC fields and 
higher laser power, photo-assisted field emission from 
higher-excitation order and photoemission processes be- 
come dominant. Together with the simulations, the im- 
portance of the electron dynamics during the electron 
emission is clarified. However, the simulations failed to 
reproduce the experimental data in terms of the emis- 
sion intensity in the higher laser power regime due to the 
neglect of the transport effects. 

A main advantage of the optical control of the field 
emission sites is that it can manipulate the electron emis- 
sion within the coherence length and time of the electrons 
inside the solid, which are roughly 100 nm and 100 fs for 
tungsten at low temperature [5|. Hence, this technique 
is maybe used to measure the coherence of the electron 
in the solids [Til. ITsT]. Moreover, the clarification of the 
detailed emission mechanism should be useful to design 
pulsed electron sources for various applications as men- 
tioned in the introduction. 



Illumination of the sharp metallic tip apex with fem- 
tosecond laser pulses creates asymmetric local fields and 
the local field distribution changes with varying laser po- 
larization. Thus, an ultrafast pulsed field-emission source 
with site selectivity on the scale of a few tens of nanome- 
ters was realized. The energy spectra of the laser-induced 
field emission reveal that the emission mechanism is dom- 
inated by photo-assisted field emission in the investigated 
regime. The emission processes and the spectral shapes 
can be tuned by the DC field and the laser power. In 
the higher DC field, photo-assisted field emission is dom- 
inant and the emission from lower-excitation order grows 
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